Geometry of large-scale low-energy excitations in the one-dimensional Ising spin glass 

with power-law interactions 
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Results are presented for the geometry of low-energy excitations in the one-dimensional Ising 
spin chain with power-law interactions, in which the model parameters are chosen to yield a finite 
spin-glass transition temperature. Both finite-temperature and ground-state studies are carried out. 
For the range of sizes studied the data cannot be fitted to any of the standard spin-glass scenarios 
without including corrections to scaling. Incorporating such corrections we find that the fractal 
dimension of the surface of the excitations, is either equal to the space dimension, consistent with 
replica symmetry breaking, or very slightly less than it. The latter case is consistent with the droplet 
and "trivial-nontrivial" pictures. 
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I. INTRODUCTION 



There have been several numerical attempts 



at 



finite temperature 



1.2.3.4.5.6.7.! 



and 



temperatur82iifliiiiiSii2iiiiiSii& to better understand 
the nature of the spin-glass state for short-range spin 
glasses. These results are generally interpreted in terms 
of the two main theories for the spin-glass phase: the 
replica symmetry-breaking (RSB) picture 
and the "droplet picture"^iiS2i22i^ (DP). RSB predicts 
that excitations involving a finite fraction of the spins 
cost only a finite energy in the thermodynamic limit, 
and that the fractal dimension of the surface of these 
excitations ds is equal to the space dimension d. This 
is in contrast to DP where a low-energy excitation 
(droplet) has an energy proportional to , where £ 
is the characteristic length scale of the droplet and 
is a positive stiffness exponent. In addition, the 
surface of these excitations is fractal with dg < d. More 
recently Krzakala and Martinji as well as Palassini and 
Young, ^ suggest an intermediate picture (called "TNT" 
for trivial-nontrivial) in which droplets have a fractal 
surface with ds < d, and their energy is finite in the 
thermodynamic limit. Which of the above pictures 
describes the spin-glass state correctly is still widely 
debated. 

The RSB and TNT pictures require two stiffness expo- 
nents for the energy of large-scale excitations. There is 
convincing numerical evidence that changing the bound- 
ary conditions (e.g., from periodic to antiperiodic), which 
induces a domain wall, costs an energy which increases as 
l*^ with 6* > 0. On the other hand, in the RSB and TNT 
pictures, the energy of droplets, created by thermal noise 
or by applying a perturbation for a fixed set of boundary 
conditions, varies as £^ with 9' = 0. By contrast, the DP 
makes the reasonable ansatz that 9' = 9 (> 0). 

In a previous publication's^ we studied the one- 
dimensional long-range Ising spin glass with power-law 
interactions. The model's advantage is that large sys- 



tem sizes can be studied, in contrast to the short- 
range spin-glass models commonly used. The results 
of Ref. showed that the stiffness exponent 9 for 
zero-temperature domain-wall excitations is positive and 
in fair agreement with analytical predictions from the 
droplet modeliSi2SiS& However, the stiffness exponent 9' 
for thermally induced droplet excitations is different and 
consistent with zero. Hence, at least for the range of 
system sizes studied, L < 512, the data of Ref. are 
consistent with both the TNT and RSB scenarios since 
they have 9> 9' = 0. 

The purpose of the present paper is to estimate dg, 
because this distinguishes between the RSB and TNT 
scenarios, since dg — d in RSB while dg < d in TNT. For 
short-range models, a droplet excitation forms a single 
connected piece, and so dg has to be zero in d = 1. How- 
ever, for long-range interactions, a droplet may consist of 
disconnected pieces,^- so a nontrivial value of dg is possi- 
ble in d = 1. We perform both finite-temperature Monte 
Carlo simulations and ground-state studies. Our results 
suggest that droplets are possibly compact in agreement 
with RSB, although the data are also consistent with a 
very small value of d — ds, which would be consistent with 
TNT. 

In Sec. m we introduce the model, observables, and 
details of the Monte Carlo technique. Results at zero 
temperature are presented in Sec. lIIII and those at finite 
temperature are presented in Sec. IIVI Our conclusions 
are summarized in Sec. Ivl 



II. MODEL AND NUMERICAL METHOD 

The Hamiltonian of the one-dimensional long-range 
Ising spin chain with power-law interactions is given by 



H — ^ JijSiSj 

{i,3) 



(1) 
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where the Ising spins Si = ±1 are evenly distributed on 
a circular ring of length L to ensure periodic boundary 
conditions. The sum is over all pairs of spins on the chain 
and the couplings are given by 
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is the straight-line distance between sites i and j. The 
random part of the interactions e^j is chosen according to 
a Gaussian distribution with zero mean and standard de- 
viation unity, and the constant c(cr) in Eq. |5J) is chosenS^ 
to give a mean- field transition temperature = 1. 

The one-dimensional long-range Ising spin chain has a 
very rich phase diagram2i2SiS& in the d-a plane. Spin- 
glass behavior is controlled by the long-range part of the 
interaction if a is sufficiently small, and by the short- 
range part if a is sufficiently large. In this work we focus 
on the long-range behavior at cr = 0.75 for whichSSi Tc > 
and the critical exponents are non-mean-field like. Using 
the exact relation^'^-^^ 6 — d — a, we expect 9 — 0.25 for 
d — I, which is in moderate agreement with numerical 
resultE^^ for domain walls induced by a change in bound- 
ary conditions at T = 0. By studying thermally induced 
droplet excitations Ref. 25 also estimated 6' « 0, consis- 
tent with RSB and TNT. 

In order to excite droplets at zero temperature we 
use the coupling-dependent ground-state perturbation 
method described elsewhere j^'^Ti^? First, we compute the 
ground-state configuration {S'{°^}. Then we perturb the 
couplings Jij by the following amount: 
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where e is a coupling constant and [• • -Jav represents a 
disorder average. The (total) energy of the unperturbed 
ground state then increases by exactly e, whereas the 
energy of any other state a will increase by the lesser 
amount of eg;, where qi is the link overlap between the 
unperturbed ground state and a state a: 



TABLE I: Parameters of the T = simulations. The ta- 
ble shows the total number of Monte Carlo steps used for 
each value of e and L. All data are computed with 10* dis- 
order realizations. The lowest temperature used to calcu- 
late the ground states with parallel tempering Monte Carlo is 
T = 0.05, the highest 1.70. We use between 10 and 23 tem- 
peratures, depending on the system size, to ensure that the 
acceptance ratios of the parallel tempering moves are larger 
than ~ 0.30. 



e 


L = 16 


L = 32 


L = 64 


L = 128 


L = 256 


L = 512 


0.50 


2 X 10-* 


4 X lO'' 


8 X lO'' 


4.0 X 10* 


3.2 X lO'' 


4.0 X 10" 


1.00 


2 X 10^ 


4 X 10'^ 


8 X 10^ 


4.0 X 10* 


3.2 X 10^ 


4.0 X 10" 


2.00 


2 X 10-' 


4 X 10^ 


8 X 10^ 


4.0 X 10* 


3.2 X lO'^ 


4.0 X 10" 



zero-temperature simulations the term "link overlap" will 
hereafter refer to the link overlap between the perturbed 
and unperturbed ground states. 

Ground states are calculated using the parallel temper- 
ing Monte Carlo method22i2£jat very low temperatures) 
as described in Refs. l25l and Isil The parameters used in 
the T — simulations are shown in Table For each 
value of L and e we compute 10^ disorder realizations. 
We find that for a — 0.75, when the model is in the long- 
range phase, the efficiency of the used algorithm to cal- 
culate ground states scales as with z = 2.9 ± 0.3. For 
the current project this translates to a total CPU time 
of 70 CPU years. Curiously, for a = 2.50, for which the 
interactions are effectively short range so frustration is 
minimal in the d — 1 model studied here, the algorithm 
performs poorly with the equilibration time varying as 
~ exp(aL), with a = 0.13 ± 0.02. It would be useful to 
understand intuitively the reasons for this. 

One quantity that we study at T = is the link over- 
lap, averaged over all samples. To see how this varies 
with size^ consider a large cluster of excited spins. This 
has a characteristic energy of order ~ , which is to be 
compared with the energy gained from the perturbation 
e(l - qi) - eL-('*-'*=). There is a distribution of cluster 
energies which we assume to have a finite weight at the 
origin, so the probability that the perturbation will cre- 
ate the excitation is ^ eL~(^ +d-ds)^ When this occurs 
1 — 9; ~ and so on average^iSLSi 



[l-qiU^eL-^'ia + bL-') , 
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where 



2^ [JjjU-v (a) MO) c-(O) 
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In previous work qi has been defined for nearest-neighbor 
models in which the sum is over nearest-neighbor pairs. 
Here we have generalized the link overlap to long-range 
models in a natural way. Because the coupling constant 
e is of order unity and not of order L only low-energy ex- 
citations can be generated. We compute the new ground 
state of the perturbed system and record the link overlap 
between the old and new ground states. In the context of 
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and we have added a correction to scaling term hL~'^. 

In RSB we have /i; = so Eq. JBJ tends to a constant 
for L ^ CX3, whereas in DP and TNT /i; > so Eq. ® 
tends to zero in this limit. 

In addition, we consider averages over only those sam- 
ples in which a large excitation is generated^^^ com- 
prising a finite fraction of spins. The criterion we take is 
l^l < 0.50. Averaging just over these samples gives^^^ 
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L-('i-ds\a + bL-") , 
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TABLE II: Parameters of the finite-T simulations. A'samp is 
the number of samples, Ns„ccp is the total number of Monte 
Carlo sweeps for each of the 2Nt replicas for a single sample, 
and Nt is the number of temperatures used in the parallel 
tempering method. 



L N samp ^ sweep Nt 

16 2.0 X 10"* 2.0 X 10^ W 

32 2.0 X 10* 4.0 X 10^ 10 

64 2.0 X 10" 8.0 X 10^ 12 

128 2.0 X lO'* 4.0 X 10" 14 

256 1.0 X 10" 2.0 X 10^ 17 

512 5.0 X 10^ 8.0 X 10'^ 24 



the prime representing the restricted average. Equa- 
tion © follows trivially from the arguments presented 
in the derivation of Eq. © with the probability factor 
+d-d,) j-gpiaced by unity. We expect that [1 — qi]'^^ 
will be independent of e. 

In order to study droplet geometries at finite temper- 
atures, we compute the distribution of the link overlap 
qi between two replicas a and (3 of the system with the 
same disorder: 

2 [JjjUv r/g(°)o(°)c(fflc(/3)\i (Q\ 

9' - ]^ pWyiK^ ^ )Jav . (9) 

Here (■ • •) represents a thermal average, and [• • -Jav rep- 
resents a disorder average. From the finite-size scaling 
arguments" used to derive Eq. © we expect that the 
variance of the distribution of the link overlap scales as 

V'Ai{qi) = L-t" {a + bL-") . (10) 

Note that in RSB ^/ = so Var(g;) tends to a constant 
for i — > oo. However, /j,; > in DP (since 9' = 9 > and 
ds < d) and in TNT (since ds < d). The bL'" term is 
a correction to scaling, which turns out to be necessary 
since the data cannot be fitted without it. 

To speed up equilibration of the finite-T simulations 
we use the parallel tempering Monte Carlo methodi22i2£ 
We test for equilibration using the criterion developed 
earlier (4 now generalized^^ for the Hamiltonian in Eq. Q • 
For all sizes, the lowest temperature used is T = 0.05, 
well below Tc ~ O.QSm^^ The highest temperature is 
1.70 which is well above the mean-field critical tem- 
perature ( = 1) and so the spins equilibrate fast 
there. We choose the spacing between the temperatures 
such that the acceptance ratios for the global moves are 
around 0.30. Parameters of the finite-T simulations are 
summarized in Table ITU 

To summarize, for L ^ cx3 all the quantities that we 
calculate, [1 — qi]a,v in Eq. JSj), [1 — qi]'^^ in Eq. ((SJ, and 
Va.r{qi) in Eq. IjlOfl tend to a non-zero constant in RSB, 
whereas they tend to zero with a power of L in TNT and 
DP. 
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FIG. 1: Zero-temperature data for [1 — g;]av as a function of 
system size L for different values of the coupling constant e. 
Note that the data only depend slightly on e, thus indicating 
only small deviations from the scaling form. The dashed lines 
correspond to a three-parameter fit to a bL~'^ as expected 
in RSB. 



TABLE III: Fits of zero-temperature data for [1 — qi]av to 
L-'-'^-'^^'>{a + bL-"), appropriate for DP/TNT, for different 
coupling constants e. The last column is P^r degree of 
freedom, where for this data with six points and four fitting 
parameters, the number of degrees of freedom (ndf) is two. 



e 


d — ds 


a 


b 


c 


xVndf 


0.50 
1.00 
2.00 


0.043(14) 
0.003(28) 
0.019(19) 


0.81(8) 

0.59(14) 

0.67(10) 


1.95(64) 

1.29(7) 

1.30(8) 


0.83(19) 

0.51(9) 

0.54(8) 


0.64 
0.31 
1.63 



III. RESULTS AT ZERO TEMPERATURE 

We first discuss the results for the constrained average 
of 1 — qi, including only samples \q\ < 0.5, since this 
yields d — ds independent of 9', see Eq. (jHJ. The results 
are shown in Fig. ^ Note that the data only depend 
slightly on e, indicating only small deviations from the 
expected scaling form. 

The results of a DP/TNT fit to [1 - qi]'^^ = 
l-id-d,)^^ + bL-") are presented in Table ITTTl while the 
corresponding RSB fits to [1 — qi]'g^ — a + bL~'^ are shown 
in Table HVI In the DP/TNT fits we find that d - ds is 
close to zero. 

Both DP/TNT and RSB fits are acceptable (xV^df ~ 
1). However, in fits to a nonlinear model, one cannot 
convert x^/ndf to a confidence limiliS even if the data 
have a normal distribution. Similarly, unlike for the case 
for fits to a linear model, the error bars do not correspond 
to a 68% confidence. We are particularly interested to get 
a confidence limit on the value of d — in the DP /TNT 
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TABLE IV: Fits of zero-temperature data for [1 — ?;]av = 
a + bL~'^, appropriate for RSB, for different values of the 
coupling constant e. The number of degrees of freedom (ndf) 
here is three. 



e 


a 


b 


c 


xVndf 


0.50 


0.580(7) 


1.35(8) 


0.53(3) 


1.74 


1.00 


0.575(6) 


1.29(5) 


0.50(2) 


0.21 


2.00 


0.569(5) 


1.28(4) 


0.474(14) 


1.36 




-0.05 

X = d-d^ 

The cumulative probability for d ■ 



0.05 



da from the fits 



FIG. 2: 

to [1 — qi]'i^v as discussed in the text. The inner pair of dashed 
horizontal lines show 68% confidence levels and the outer pair 
show 95.5% confidence levels. 
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FIG. 3: Zero-temperature data for [1 — g;]av as a function 
of system size L for different values of the coupling constant 
e. The dashed lines represent fits according to [1 — g;]av = 
a + bL-" (RSB). 



TABLE V: Fits of the zero-temperature data for [1 — g;]av = 
L"^' (a + 6L"''), which assume the DP/TNT picture, for dif- 
ferent coupling constants e. 



e 


Mi 


a 


b 


c 


xVndf 


0.50 
1.00 
2.00 


0.065(62) 
-0.15(18) 
0.018(7) 


0.28(1) 

0.07(12) 

0.44(25) 


1.0(8) 

1.07(7) 

1.25(8) 


0.81(46) 

0.51(6) 

0.49(13) 


0.12 
0.68 
0.64 



fits. We do this by computing a-s a function of — d^, 
minimizing with respect to the other parameters (a, fo, 
and c). The probabihty of the fit P is proportional to 
exp(— which we numerically integrate to get the 
cumulative probability for x ~ d ^ ds'. 

QM^fpWW^ (11) 

The results are shown in Fig. |21 The data for e = 0.5 
and 2.0 constrain d — dg to zero or a small positive value. 
The data for e = 1.0 constrain d ~ dg less and allow a 
range of negative values which are unphysical. At a 68% 
confidence level the data are consistent with 

0<d-4<0.05, (12) 

apart from the e = 0.5 data which would exclude zero 
at the 68% level but, from Fig. [3 are consistent with it 
at the 86% level. We take Eq. H12|l to be our estimate 
for d — ds- It is consistent with the RSB prediction of 
zero and also consistent with a small non-zero value in 
the DP/TNT scenarios. 

Data for the average oil — qi over all samples are shown 
in Fig. O along with RSB fits. The data show curvature 



indicative that corrections to scaling have to be included. 
The DP/TNT fits, Eq. ®, are shown in Table |V| while 
the RSB fits (which fix /i/ to zero) are shown in Table 
IVII Both fits have acceptable . 

The cumulative probabilities shown in Fig.^give a lot 
of weight to unphysical negative values of /ij, especially 
for e = 1. For all values of e the weight is small for 
greater than about 0.10 so we conclude that 

0<^i<0.10. (13) 

We should, perhaps, be cautious about this statement 
in view of the large weight at negative values of /i; in 
Fig.^ However, Eq. ifT^ is consistent with Eq. ifT^ and 
the resuh of Ref. ^ that 0' ~ 0. 



TABLE VI: RSB fits of zero-temperature data for [1 — gjjav = 
a + bL^'^ for different values of the coupling constant e. 



e 


a 


b 


c 


xVndf 


0.50 


0.168(8) 


0.73(11) 


0.56(7) 


0.25 


1.00 


0.280(8) 


1.15(11) 


0.57(5) 


0.90 


2.00 


0.378(10) 


1.25(6) 


0.46(3) 


0.45 



X 




FIG. 4: The cumulative probability for d ~ ds from the fits 
to [1 — gjjav as discussed in the text. The inner pair of dashed 
horizontal lines show 68% confidence levels and the outer pair 
show 95.5% confidence levels. 
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FIG. 5: Log-log plot of finite-T data for the variance of the 
link overlap Var(g; ) as a function of system size L for several 
low temperatures. In all three cases we see strong curvature in 
the data suggesting corrections to scaling. The dashed lines 
represent fits according to a + hL~'^ (RSB) with the fitting 
parameters shown in Table IVnTI 



IV. RESULTS AT FINITE TEMPERATURE 

In this section we study the model at temperatures well 
belowS^ Tc « 0.63. Figure shows data for the variance 
of the link overlap for several low temperatures. The 
data show strong curvature indicative that a simple fit of 
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Ml 


0.05 


-0.21 


0.10 


0.10 


0.16 


0.16 


0.23 


0.13 




0.05 



0.1 
Ml 



0.15 
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FIG. 6; as a function of ^i, optimized with respect to the 
other parameters (a, 6, and c) in Eq. IIOII . for the variance of 
the link overlap. The arrows mark the minima. 



TABLE VIL DP/TNT fits of Var(gi) to L"''' (a + bL-") for 
different temperatures. 



T 


Mi 


a 


6 


c 


xVndf 


0.05 


-0.21(54) 


0.002(11) 


0.073(8) 


0.52(35) 


0.05 


0.10 


0.10(11) 


0.047(42) 


0.16(3) 


0.55(24) 


0.52 


0.16 


0.16(6) 


0.079(39) 


0.29(6) 


0.60(17) 


0.83 


0.23 


0.13(8) 


0.050(34) 


0.41(2) 


0.52(2) 


1.18 



the form aL ^' is improbable and corrections to scaling 
must be included. 

Fits to Eq. 10 (DP/TNT picture) are shown in Table 
IVIII and fits to the RSB picture (in which /i; is fixed to 
be zero) are shown in Table IVlTll The quality of the fits 
is acceptable. 

However, the DP/TNT fit for T = 0.05 gives an un- 
physical negative value for fii with a very small amplitude 
a. To clarify this situation, we plot, in Fig. as a 
function of optimizing with respect to the other pa- 
rameters (a, 5, and c). For T = 0.05, is quite small 
out to very large negative values of /z/ (not shown) and 
increases rapidly for ^/ greater than about 0.12. Since 
physically /i/ cannot be negative, the only conclusion we 
can deduce from the T — 0.05 data is that m lies between 



TABLE VIII: 
temperatures. 



RSB fits of Var(g;) to a + bL-" for different 



T 


a 


b 


c 


xVndf 


0.05 


0.015(2) 


0.073(10) 


0.47(7) 


0.10 


0.10 


0.020(2) 


0.155(11) 


0.47(3) 


0.47 


0.16 


0.021(1) 


0.261(11) 


0.50(2) 


1.32 


0.23 


0.017(1) 


0.391(9) 


0.55(1) 


1.35 
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zero and about 0.12, consistent with the result from the 
T = data in Eq. ||T3J|. The data for for T = 0.10 in 
Fig. El has a minimum at /i/ — 0.10 but it is shallow and 
fii — has only a slightly greater value. The T = 0.10 
data are therefore also consistent with Eq. H13|l . The data 
at higher temperatures, T = 0.16 and 0.23 have a deeper 
minimum at nonzero x^i suggesting that /i; = is some- 
what unlikely, but experience from short-range systems^ 
suggests that estimates of m at finite T are effective ex- 
ponents which need to be extrapolated to T = to get 
close to the asymptotic value. Hence we do not feel that 
the results at T = 0.16 and 0.23 rule out /i/ = 0. 

Overall, the finite-T data are consistent with /i/ in the 
range given by Eq. H13|l which came from the T = data, 
and do not constrain /i/ any further. 

V. CONCLUSIONS 

We have studied the geometry of the large-scale, low- 
energy excitations in a one-dimensional Ising spin glass 
where the interactions fall off as r^"' with a — 0.75, both 
at r = and at temperatures well below the spin-glass 
transition temperature. We find that the fractal dimen- 
sion of the surface of these excitations, ds, lies in the 
range < c? — < 0.05. This is consistent with the 
RSB picture (d — = 0). It is also consistent with the 
DP/TNT picture {d - ds > 0) but with a small value 



of d — dg. Substantial corrections to scaling had to be 
incorporated into all the fits. 

We have also estimated the exponent /i/ — 6'+2{d—ds), 
where 9' characterizes the dependence of the energy of 
droplet excitations on their length scale. We find it to be 
in the range < /i; < 0.10, which is consistent with the 
value for d — d^ in Eq. H12|l and our earlier result^^ that 
6*' ~ 0. Note that this result for 6*', if also valid in the 
thermodynamic limit, is inconsistent with the DP. 

By studying a one-dimensional model, we have been 
able to study a much larger range of sizes, 16 < L < 512, 
than is generally possible in spin glasses. However, in the 
absence of a good understanding of corrections to scaling 
in spin glasses, we still cannot rule out the possibility 
that different behavior may occur in the thermodynamic 
limit. 
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